Thermodynamic transitions in inhomogeneous cuprate superconductors 
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Scanning tunneling spectroscopy studies on Bi2Sr2CaCu2 08+« suggest the presence of electronic 
inhomogeneity with a large spatial variation in gap size. Andersen et al have modelled this variation 
by assuming a spatially-varying pairing interaction. We show that their calculated specific heat is 
incompatible with the experimental data which exhibit narrow transitions. This calls into question 
the now-common assumption of gap and pairing inhomogeneity 

PACS numbers: 74.25.Bt, 74.40.4-k, 74.81.-g, 74.72.-h 



In a recent paper Andersen and coworkers^ at- 
tempt to model the spectroscopic inhomogeneity inferred 
from scanning tunneling spectroscopy (STS) studies on 
BizSrzCaCuaOg+a (Bi-2212)2.3. These STS studies sug- 
gest a ±25% variation in both the gap magnitude and 
the local density of states (LDOS), on a length scale of 
the order of one or two coherence lengths, <^o- The loca- 
tions with large spectral gaps lack well-developed coher- 
ence peaks and, based on the absence of Ni resonances 
there, they are suggested to be non-superconducting-. 
These conclusions have become highly influential. 

To model these effects Andersen et aZi adopt a 
spatially-inhomogeneous pairing interaction in a d-wave 
BCS pairing model and solve self-consistently for the lo- 
cal gap magnitude in order to map the spatially varying 
interaction onto the observed variation in gap magnitude. 
They then compute the specific heat and show that the 
breadth of the anomaly is similar to the spread in local 
gap magnitude. This breadth, they claim, is "compara- 
ble to experiments on this material" , and conclude that 
"substantial nanoscale electronic inhomogeneity is char- 
acteristic of the bulk BSCCO system" . 

However our data^ (see Fig. 1), which they use for 
their comparison, points to the very opposite conclusion. 
They erroneously equate the extended fluctuation region 
above Tc with broadening of their mean-field (MF) tran- 
sition over a 40K range. Strong fluctuations are an in- 
trinsic property of this highly-anisotropic material^ but 
are not included in their calculations. In fact it is the 
narrow region of strong negative curvature close to each 
peak (see arrows in Fig. I) that reveals the true extent of 
extrinsic broadening. Such a strong T-dependence over 
a very restricted T-range would not be possible in a ma- 
terial with a broad distribution of Tc values. 

Moreover, the locations where their gap is a maximum 
correspond to the maximal local pairing interaction i.e. 
where SC is strongest. The STS studies'^ show the op- 
posite: SC is weakest and perhaps absent at the points 
where the supposed gap is maximal. Their model also im- 
plies that there are local regions where SC persists well 
above as much as 30K for Bi-2212. Such regions are 
not observed in STSSii. 



Here we consider the experimental specific heat data in 
more detail and show that the transitions are not strongly 
broadened as suggested by Andersen et a/i, thus poten- 
tially reversing their inference of inhomogeneous pairing 
in the bulk material. This concurs with several other 
studies which imply the absence of gross inhomogeneity 
in cuprate superconductors^iS. 

The specific heat near Tc consists of a MF step at 
Tc and a (nearly symmetrical) fluctuation contribution 
above and below T^. More generally, Tc may be broad- 
ened out into a distribution of Tc values. The separate 
contributions of fluctuations and transition broadening 
may seem similar well away from the mean Tc, but nearby 
they are quite distinctive and easily separated. Andersen 
et ali fail to recognize this distinction. 

Let us consider first fluctuations where there is a 
sharply defined Tc. The fluctuation specific heat should 
diverge at Tc but is cut-off due to the inhomogeneity 
length scale. For Bi-2212 we have previously analyzed the 
fluctuation contribution* and deduced transition half- 
widths as small as ATc/Tc ^ 0.014, consistent with an 
inhomogeneity length scale as large as 16^0i much greater 
than ^0- The cutoff is reflected in the narrow region of 
negative curvature between the inflexion points in the 
specific heat coefficient, 7(T), just above and below Tc 
as shown by the arrows in Fig. 1(a). 

In Fig. 1(b) we show the derivative d^/dT from our 
data and compare it with that from Andersen's model 
calculation. The data curves correspond to doping states 
of p = 0.16, 0.19, 0.20 and 0.215. The inflexion points 
are located at the maxima and minima below and above 
Tc, and between them d^jdT changes sign. For -p ~ 0.19 
the inflexion points are just 3.3K apart. For the model 
calculation they are up to 40K apart (dashed curve and 
arrows), just as would be expected for a ±25% spread 
of gap values. This directly illustrates that we observe 
rather sharp thermodynamic transitions in Bi-2212. 

To show that the transition is indeed narrow we plot 
in Fig. 2 the field dependence of 7 for a Bi-2212 sam- 
ple with doping p « 0.21. Firstly it is clear from this 
plot that Tcili = 0) is close to the peak (as expected if 
the MF step is small) . The narrow peak is progressively 
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FIG. 1: (a) The specific heat coefficient, 7, for Bi-2212 with 
p= 0.16, 0.19, 0.20 and 0.215, respectively, (b) The derivative 
d^/dT. The dashed curve is dj/dT from the Andersen et al. 
calculation. Arrows indicate inflexion points in 7(T). 



FIG. 2: (a) The field dependence of the specific heat coef- 
ficient 7-7(ref) for Bi-2212. The field is shown in units of 
Tesla. (b) the derivative, ^ x (7 — 7(ref)) showing the field- 
broadening of the inflexion points. 



suppressed and broadened by the field, with a marked ef- 
fect even for fields as low as 0.3T. The vortex separation 
^J(t>Q/H ~ 45 nm/i/i/ (Tesla), which acts as an inho- 
mogeneity length scale, is very large at such low fields 
and the sensitivity of the transition to a field as low as 
0.3T supports our conclusion that the order parameter 
is rather homogeneous. Calculated transitions based on 
an inhomogeneity length scale of 1.6 - 2.5 nm would be 
totally insensitive to such low fields. 

To gain more insight into the shape of the anomaly we 
have simulated the effects of broadening by integrating 
sharp specific heat transitions over gaussian distributions 
of Tc- The data are compared with similar plots for some 
of our Bi-2212 data. To cover a wide range of situations 
we have considered (a) transitions with a simple MF step, 
(b) pure Int fluctuations with no MF step, and (c) ad- 
mixtures of the two. In all following cases C corresponds 
to the difference C — C normal- 

For the simulated MF specific heat function we adopt 
a form which conserves entropy at and has a step 
AC„/ = 2 at To. 

Cmf{x) = x(3x2 - 1) for X = T/To < 1, , . 
Cmf ix) =0 for a; > 1. ^ ' 

For the fluctuation specific heat we assume 

Cfiuc{x) = ln(l/ I a; - 1 I) for all x (2) 

This has a symmetrical divergence at T^. For the dis- 
tribution of Tc values we assume a Gaussian function 



centered on Too with half-width w = 0.001, 0.02, 0.05, 
0.1, 0.15 and 0.2: 

The specific heat resulting from this Tc distribution is 

poo 

C{T) = / C{x)P{y)dy (4) 

Results. For a range of widths w, plots of C and dC/dT 
versus t = [T/Tco) — 1 are shown in Fig. 3 for Cmf and 
in Fig. 4 for C/j„c • Admixtures Cmix = Cmf + mCfiuc 
with TO = 0, 0.2, 0.4, 0.6, 0.8 and 1 are plotted versus t 
in Fig. 5(a) and (b) for w = 0.02. 

It is important, firstly, to note that in spite of the very 
different T-depcndencics of the unbroadcned anomalies, 
the region over which broadening effects arc important 
is very similar for Cmf, Cfiuc and Cmix. This shows 
that, though the Andersen model does not include fluc- 
tuations, the addition of fluctuations would not alleviate 
the disparity between their model and the experimental 
data. In particular, inclusion of fluctuations will not re- 
sult in a more narrow negative-curvature region around 
Tc as seen in the data. We also find that the relative posi- 
tions of key features of the T-derivatives of the broadened 
curves are insensitive to the detailed T-dependencies as- 
sumed for Cmf and Cfiuc. From these plots we find the 
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FIG. 3: (a) The MF specific-heat anomaly modeUed using 
eqn. (1) with a normal distribution of Tc values with half 
width, w = 0.001, 0.02, 0.05, 0.10, 0.15 and 0.20. (b) The 
derivative dCm//dr. Inset: 5t (defined in text) versus w. 



FIG. 4; (a) The fluctuation specific-heat anomaly modelled 
using eqn. (2) with a normal distribution of Tc values with 
half width, w = 0.001, 0.02, 0.05, 0.10, 0.15 and 0.20. (b) The 
derivative ACfiuc/AT. Inset: 5t (defined in text) versus w. 



locations of key features and determine their relation to 
the half- width w of the distribution of values. 

(i) For the MF anomaly, shown in Fig. 3, the most 
useful features are the positions of the points of maximum 
negative and positive curvature (PC/dT^ at i_ and 
respectively, and the difference 5t — t+ — t-. In the inset 
to Fig. 3(b) we show 5t plotted as a function of the 
half-width w. Over most of the range of w, 5t k, 2.Qw. 

(ii) For the symmetric fluctuation anomaly, shown in 
Fig. 4, the most useful features are the positions of the 
points of maximum positive and negative slope dC / dT at 
t- and t+ respectively, and the difference St = — t^. 
In the inset to Fig. 4(b) we show St plotted as a function 
of w. Over the entire range, St « 2.63w. 

(iii) For the admixture of a MF anomaly and fluctua- 
tions. Fig. 5 shows plots of Cmix = C'mf + rnCfiuc for 
w=0.02 with m = 0, 0.2, 0.4, 0.6, 0.8 and 1. These are 
useful for comparison with typical cuprate speciflc heat 
data and cover the range from asymmetric (m=0) to al- 
most symmetric (m=l) anomalies. A value of to « 0.8 
is appropriate for weakly-overdoped Bi-2212. From Fig. 
5(b) we obtain the positions of the points of maximum 
positive and negative slope at t^ and respectively, 
and the difference St = t+ — t^. For larger values of to, 
typical of our Bi-2212 samples, we find St « 2.7w. For 
TO=oo (pure fluctuations) we found, above, St « 2.63w. 

Now we compare this model data with measured data 
for Bi-2212. Plots of dC/dT versus r = T/Tp - 1 are 
shown in Fig. 5(c) for eight values of p (ranging from 



well underdoped to well overdoped). Comparison with 
Fig. 5(b) shows that Tc is only 1-2% above Tp, so sepa- 
rations of peaks in t and in r are almost identical. When 
underdoped the anomaly is pure fluctuations with no MF 
step while the increase in MF step is evident for overdop- 
ing from the increasingly different relative magnitudes of 
the positive and negative peaks in dC/dT. 

We estimate the transition width, w, from the separa- 
tion of the positive and negative peaks in dC/dT. For 
all samples with p > 0.169 the separation is St « 0.05 - 
0.06. Taking St sa 2.7w for to f^l (see Fig. 5(b)) gives 
w « 0.019 - 0.021 for the half width of the distribution 
of Tc values. For p=0.162 we have St w 0.073 or w w 
0.028, and for ]3=0.138 we have St « 0.14 or w « 0.056. 
All these values are in good agreement with the values of 
the half-width At shown in Fig. 4 of our previous work^. 

As a check on our previous method* for estimating 
these half-widths. At, we show in Fig. 6 plots of Cmix 
versus logiQ(i) and log^ol^*), respectively, for w = 0.02 
and to = to 1, where t* = Vt^ + A?. Our normal 
procedure is to choose a value of At that just averts the 
negative curvature close to Tc seen in plots of Cmix versus 
logiQ{t), and it is evident from the solid curves that the 
choice At — w achieves this result. This confirms that 
estimates of the broadening At from plots of AC ver- 
sus logiQ(t*) give reliable values comparable to the true 
half-width w. The plots in Fig. 6 also show that this 
procedure provides a reliable estimate of the MF step 
ACmf ~ 1.5 even in the presence of large fluctuations. 
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Turning to the MF calculation of the specific-heat 
anomaly by Andersen et al^, and recalling that their 
model does not include fluctuations, we compare their 
results (shown in Fig. 4(b) of their paper) with our plots 
for a fluctuation-free broadened MF step, shown here in 
Fig. 3. From the location of the temperatures T_ and 
T+ of their maximum negative and positive curvature, 
we obtain 6t = {T+ - T_)/Tav « 0.146, 0.204 and 0.232 
for their curves for 6g/t = 1.0, 1.5 and 2.0, respectively. 
Taking 6t = 2w from Fig. 3(b) we find w « 0.073, 0.102 
and 0.116, respectively, as the half-widths for the three 
curves. This agrees with the « 20% breadth of the tran- 
sition that they quoted for dg/t = 1.5, the value which 
best accounts for the gap-maps. This may be compared 
with the similar analysis of our Bi-2212 data from which 
we obtained w « 0.02 over most of the doping range, 
a factor of five, or more, lower than the Andersen esti- 
mate based on gap-maps. It surely cannot be claimed 
that the spectroscopic and thermodynamic data are con- 
sistent. Indeed, given our evidence that the spread of Tc 
values is rather narrow, their calculations show plainly 
that the inhomogeneity in the gap-maps cannot result 
from gross pairing energy disorder in the bulk. 
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FIG. 5: (a) The admixture specific-heat anomaly Cmix ~ 
Cmf + rnCfluc with half width, w — 0.02 and various mixing 
ratios, m, plotted versus t = (T/Tco) — 1- (b) The derivative 
dC'mix/dT vs t. The zero crossing occurs at Tp and the curves 
intersect at Tco. (c) dC/dT vs r for Bi-2212 at diS'erent doping 
levels (annotated) showing a crossover from pure fluctuations 
for p < 0.162 to an admixture with a MF step for p > 0.162. 
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FIG. 6: The admixture specific-heat anomaly Cmix = Cmf + 
mCfiuc with half width, w = 0.02, plotted against logj^Q(f) 
(dashed curves) and against logjQ(t*) (solid curves) with t* — 
\/t^ + At^ and setting At — w. 



To conclude, we have shown that the inference of pair- 
ing inhomogeneity from STS gap maps, and the resultant 
transition broadening, is inconsistent with the specific 
heat data which exhibit sharp features with transition 
widths of the order of 3K in Bi-2212. It is not possi- 
ble with any broad spread of SC gaps to have strong T- 
dependences over a narrow T-range. This concurs with 
a recent analysis of quasi-particle scattering seen in spa- 
tial modulations of STS data. The Fourier transform 
of these patterns reveals spots, corresponding to scatter- 
ing 5- vectors, that are far too narrow for the presumed 
broad distribution of SC gaps^. The inference of large- 
scale gap and pairing inhomogeneity at (7r,0) does not 
appear to be supported by the wider thermodynamio^, 
NMM, ARPESia and tunneling2*ii data. 
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